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Abstract 

We present a self-consistent approach to the modehng of dense plasma mixtures in local 
thermodynamic equilibrium. In each electron configuration the nucleus is totally screened by 
electrons in a Wigner-Seitz sphere (ion-sphere model) . Bound and free electrons are treated 
quantum-mechanically. The assumption that all species should have the same electronic 
environment leads to the equality of the electronic pressure for all ions of all elements having 
therefore different cell volumes. The variation of the average atomic radii of the different 
elements with respect to temperature is investigated, and the procedure is applied to the 
determination of pressure in the Sun center. 

1 Introduction 

Theoretical studies of electronic structure of disordered systems composed by ions and electrons 
are of great interest for understanding of dense-plasma radiative transfer and equation of state 
(EOS) [UEII^. In astrophysics, EOS of plasma mixtures is essential for stellar pulsation calcu- 
lations, via pressure modes (requiring sound velocity given by the second derivative of the free 
energy with respect to matter density) or via gravity modes, highly sensitive to the temperature 
at the star center. Moreover, EOS has an impact on the photo-absorption cross-sections, which 
determine temperature distribution inside the star. Such EOS studies concern atoms in extreme 
conditions of density and temperature, which implies existence of multicharged ions close to local 
thermodynamic equilibrium (LTE). An ion configuration is defined by sub-shells (orbitals) occu- 
pied by an integer number of bound electrons. For example c = (ls)^(2s)^(2p)^(3s)^(3p)^(3(i)^ is 
a configuration consisting of six sub-shells containing respectively 2, 2, 5, 1, 3 and 2 electron(s). 
In our model, the nucleus is totally screened by bound and free electrons in a sphere (ion-sphere 
model) which is usually named, refering to solid-state physics, a Wigner-Seitz (WS) sphere. 
Throughout the paper, the term "element" means "element of the periodic table" (chemical 
element): H, He, Li, Be, B, C, • • • . 

In the usual so-called ideal-gas mixing rule, the electronic structure of each element is cal- 
culated separately and a mixture is obtained by simply including all the individual elements 
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Figure 1: Variation of atomic radii of carbon and oxygen with respect to the temperature 
and comparison with the results of Yuan [6]. The non-monotonicity illustrates the competition 
between electronic structures of the two atoms. 

at fixed temperature with their abundances. There are also other ways to treat mixtures, for 
instance in the framework of orbital- free molecular dynamics [1] . 

The present model relies on the idea that at the LTE conditions, all ionic species shall be 
considered at the same plasma electronic environment (free electron " bath" ) . In section [21 we 
show that, provided that ions of different charges have different volumes, the electronic pressure, 
calculated together with ionic volumes in a self-consistent way, must be equal for each ion of 
each element. Comparisons with chemical-picture [5] and average-atom mixture models [6] in 
the case of a CO2 plasma are presented in section [31 In the same section, values of pressure 
at the Sun center are compared with the results of OPAL physical-picture model [3 [HI [9] and 
Quantum Langevin Molecular Dynamics (QLMD) simulations |10j . 

2 Self-consistent-field calculation of atomic structure and equal- 
ity of electronic boundary pressure 

Atoms in a plasma can be idealized by an average atom confined in a Wigner-Seitz (WS) sphere, 
which radius r^^ is related to matter density. The electronic pressure Pg [12] consists of three 
contributions, Pe = Pb + Pj + Pxc, where the bound- and free- electron pressures, respectively 
Pf) and Pf, are evaluated using the stress-tensor formula 
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Us representing the radial part of the wavefunction ipg (corresponding to energy e^) multiplied 
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by r, Eq the rest-mass energy of the electron and i the orbital angular momentum. / denotes 
the Fermi-Dirac distribution: 
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/i being the chemical potential. For bound states, s = n,i and J = Yl'^=i J21=o ' ^^'^ 
for free states, s = e,i and = degJ^T^o- "^^^ choice of the boundary conditions plays 
a major role in the pressure value. This comes from the fact that the energy of an orbital 
depends on the value of the corresponding wavefunction at the boundary of the WS sphere. 
In our model, the bound wavefunction behaves like a decreasing exponential at the boundary. 
Pxc is the exchange-correlation pressure evaluated in the local density approximation using 
the formulas of Ichimaru et al. [13]. The adiabatic approximation is used to separate the 
thermodynamic functions into electronic and ionic components. The total pressure can be 
written P = Pi + P(, = Pi + Ph + Pf + Pxc-, where Pi is the ionic pressure, calculated according 
to [14] (ideal gas with One-Component-Plasma corrections). The self-consistent calculation of a 
configuration provides the potential, the electron populations of the orbitals and the free-electron 
chemical potential /i, resulting from the neutrality of the plasma. 

In the plasma, all ions should have the same electronic environment. All configurations of 
the same charge are grouped together and constitute an ion. Normalization of configuration 
probabilities Wc = 1 associated to Lagrange multiplier B) and preservation of total volume 
WcVc = y associated to Lagrange multiplier P) must be ensured. The constrainted grand 
potential is: 
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and the variation with respect to volume Vc and probability Wc give respectively 
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Equalities of Eq. @ mean [16] that all configurations have the same pressure and 
that it is necessary to consider enthalpy instead of the free energy. In practical, all pressures 
are calculated in a self-consistent way and equalized by a multi-dimensional Newton-Raphson 
method. The average value of thermodynamic quantity X{e,i) for ion i{e) of element e and its 
probability Wj(e) are given by: 
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where Wc and X{c) represent respectively the probability and the value of quantity X for 
configuration c. The average value {X) of a quantity X thus reads: 
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where E is the number of elements in the plasma and the number of charge states 
of element e. In the present work, pressure equality means the equality of the stress-tensor 
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T 

(eV) 


Total pressure 
Potekhin and Chabrier 
(Mbar) 


Total pressure 
This work 
(Mbar) 


50 


2.22 


1.86 


100 


5.06 


4.46 


500 


27.16 


27.17 



Table 1: Total pressure for a CO2 plasma at p(=0.1 g/cm^. Comparison with values obtained 
from the chemical-picture model of Potekhin and Chabrier [5]. 



Element 


Z* 


P 


Pe (Mbar) 


Element 


Z* 


P 


Pe (Mbar) 


H 


1 


105.81 


2.10 10^ 


P 


13.84 


248.12 


8.50 10"^ 


He 


2 


211.74 


1.03 10^ 


S 


14.62 


243.73 


8.67 10^ 


C 


5.69 


228.02 


9.64 10"^ 


CI 


14.63 


269.70 


8.29 10^ 


N 


6.62 


229.25 


9.63 10"* 


Ar 


15.29 


291.31 


7.78 10^ 





7.57 


229.68 


9.61 lO'' 


Ca 


16.86 


265.98 


8.47 10^ 


Ne 


9.47 


233.11 


9.47 10* 


Ti 


18.52 


290.34 


7.77 10^ 


Na 


10.41 


242.45 


9.11 lO'' 


Cr 


20.23 


289.58 


7.77 10^ 


Mg 


11.33 


236.13 


9.35 lO'' 


Mn 


21.09 


293.99 


7.62 10^ 


Al 


12.23 


243.33 


9.01 lO'' 


Fe 


21.95 


287.63 


7.76 10^ 


Si 


13.13 


236.62 


9.20 lO'' 


Ni 


23.21 


281.42 


7.73 10^ 



Table 2: Partial ionizations Z* and densities p (in g/cvo?) obtained from our mixture model 
with Grevesse composition [19j. Pe represents the pressure of each element before applying the 
mixture procedure. 



pressure at the boundary of all ionic spheres. Similarly to the single-element case \Vb\ I16|. it is 
convenient to write the system of equations in the following form: -Pi(i) = -f2(i)) -^2(1) = -^3(1)' 

■■■) Pnx{1) = -^1(2)' -^1(2) = -^2(2); ^2(2) = -^3(2)5 '■■) -^^2(2) = -^1(3)' ^liE) = P2(E)^ 

P/v^_i(£;) = Pne(E)-> where Pj(e) represents pressure at the boundary of ion i of element e. The 
partial density of this ion is Pi(^e) a-nd the index i{e) varies from 1 to A'^g- The closure relation of 
the system is the matter conservation law 

where pt is the matter density of the plasma and Me and Xe represent respectively the atomic 
mass and the mass fraction of element e. The Jacobian matrix can be inverted analytically since 
only diagonal-, first super-diagonal- and last-line coefficients are non zero. 



4 



Pt (g/cm^) 


T(eV) 


EOS model 


P (Mbar) 


141.25 


989.45 


OPAL 


1.6162 10^ 






QLMD 


1.6090 10^ 






This work 


1.6217 10^ 


141.25 


1189.58 


OPAL 


1.9337 10^ 






QLMD 


1.9122 10^ 






This work 


1.9403 10^ 


152.70 


1352.64 


QLMD'* 


2.3525 10^ 






This work'* 


2.3584 10^ 






QLMD" 


2.2575 10^ 






This work** 


2.2629 10^ 






QLMD'^ 


2.2328 10^ 






This work'^ 


2.2532 10^ 



Table 3: Comparison of central pressures of the Sun at typical temperatures (T) and densities 
(p) between our approach, the OPAL equation of state [Zl El E] and the Quantum Langevin 
Molecular Dynamics [10]. (") corresponds to the mass fractions of H, He and C: X=0.3387, 
Y=0.6613 and Z=0, ('') corresponds to the mass fractions of H, He and C: X=0.3125, Y=0.6406 
and Z=0.0469 and {^) corresponds to the mass fractions of H, He and O: X=0.3077, Y=0.6308 
and Z=0.0615. 




Figure 2: Schematic picture of the statistical mixture model. 
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3 Comparisons with other approaches and appUcation to the 
Sun center 



Figure [T] shows, in the case of a CO2 plasma at a density of 0.1 g/cm^ [6], that the atomic 
radii of carbon and oxygen exhibit a non-monotonic variation with the temperature. This is 
the signature of the shell structure, absent from Thomas- Fermi model [T7]. The competition 
between the ionization of the two atoms can be seen very clearly. One can find that there are 
two main structures around 10 and 50 eV, respectively. The first increase of the carbon radius 
around 10 eV is due to the ionization of the L shell (2s and 2p orbitals). For oxygen, this occurs 
later due to the fact that the binding energies are higher than for carbon. The second increase 
of radii, around 50 eV, is due do the ionization of the K shell (Is orbital) which is ionized more 
rapidly for carbon than for oxygen. The radii obtained from the average-atom mixture model 
of Yuan [6] differ from ours mainly at low temperatures. Total pressures are compared with the 
values of Potekhin and Chabrier [5] for three different temperatures in table [TJ 

The solar neutrino flux predicted by a solar model depends on the Rosseland opacity of the 
stellar material in the Sun center [181. Opacity influences the details of nuclear reactions and 
therefore the solar luminosity. Heliosismology provides a large array of accurately measured 
(0.01 %) p-mode oscillation frequencies, related to the equation of state of the solar plasma. 
Temperature is taken to be Te=l.3621 keV and density pt=157.02 g/cm^ and the abundances of 
different elements are taken from the well-known Grevesse mixture [19]. The electron pressure 
predicted by our model is 1.4 Mbar and the total pressure 2.43 Mbar (if ions are treated as a pure 
ideal gas) and 2.46 Mbar (if OOP corrections are included). The resulting partial ionizations 
and densities are displayed in table [21 and comparisons with the OPAL values [3 [HI [9] and 
with the QLMD results of Dai et al. in table [3l The OPAL equation of state belongs to the 
"physical-picture" models, and relies on an activity expansion. The bound states are described 
through the Planck-Larkin partition function, and the systematic expansion from one-body to 
many-body terms allows one to take into account the continuum contribution (scattering states) 
|20j . Solar mixture can be a good test for our models: for instance very small quantities of 
Fe and Ni provide a large contribution to the opacity, which influences the details of nuclear 
reactions (and therefore the solar luminosity) as well as the predicted neutrino flux |18j . 

4 Conclusion 

A new approach for the thermodynamics of plasma mixtures has been presented. All electrons 
are treated quantum-mechanically and the shell structure has an important impact on the partial 
densities of different elements. The electronic pressure is imposed to be equal at the boundary 
of each ion sphere. Thus, the partial density of each ions charge state is determined in a 
self-consistent way. Future work includes further comparisons with other approaches, such as 
chemical-picture models [51[2I] and improvement of thermodynamic consistency (see Refs. P3[2])- 
We plan also to study other thermodynamic quantities, such as internal energy, sound speed 
or gravitational acceleration, and to apply the present model to the calculation of electrical 
resistivity [22] of plasma mixtures. 
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